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Abstract: We calculate the quark-anti-quark contribution to the next-to-leading log- 
arithmic corrections to the BFKL kernel, retaining the dependence on the momenta of 
the produced particles. This allows us to study the details of the NLL corrections. We 
demonstrate that the standard calculation of the NLL corrections to the scattering of two 
off-shell gluons includes contributions from energies far above that which is probed at LL. 
This explicitly violates energy and momentum conservation in the evolution and could be 
a source of the reported large NLL corrections. The presented calculation is a step towards 
combining energy and momentum conservation with full NLL accuracy in the evolution. 
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1. Introduction 

The BFKL formahsm is currently used both in the description of the small-a; evolution 
of the pdfs, and to obtain the leading logarithmic approximation (LLA) to the scattering 
matrix element for e.g. pure multi-jet and forward PF+multi-jet production at hadron 
colliders. Even though the transverse scales involved in the two problems arc very different, 
both applications rely on the same BFKL equation in the description of the evolution 
of an off-shell i-channel gluon by the emission of several gluons, and quark-anti-quark 
pairs at next-to- leading logarithmic accuracy (NLL A). Since the gluons emitted from the 
BFKL evolution are relatively hard when the formalism is applied to the approximation 
of the hard matrix element, it proves absolutely crucial to supplement the solution of the 
BFKL equation with energy and momentum conservation. While this is possible using 
the LL BFKL kernel, the derivation of the NLL BFKL kernel already contains phase space 
integrals over the invariant mass of two emitted gluons and a quark-anti-quark pair. In the 
standard analysis [1,2], the NLL corrections are found to be size-able, and it therefore seems 
imperative to include them in the calculation of any quantity which is to be confronted 
with reality. 

In this paper we will recalculate the quark-contribution to the NLL BFKL kernel in 
a form that will allow us to combine the evolution to full NLL accuracy with energy and 
momentum conservation [3] . We will use the obtained result to study the make-up of the 
NLL corrections. Specifically, we will show that in the standard analysis, due to the lack of 
energy and momentum conservation, the NLL corrections to a Reggeon-Reggeon scattering 
at a given energy receives contributions from significantly larger energies. This could be a 
significant and unintended source of the observed large NLL corrections. 

In the interest of clarity and conciseness in our argumentation, we have chosen to 
include all intermediate results in the derivation of the NLL BFKL kernel. 
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2. Exclusive States to Next-to-Leading Logarithmic Accuracy 



We will begin this section by recalling the results of the calculation of the NLL corrections 
to the BFKL kernel. In the following sections we will recalculate the quark-anti-quark 
contribution using a phase space slicing method, which allows us to study directly the 
energy dependence of the corrections to the Reggeon-Reggeon collision. 

The BFKL equation governing the evolution in rapidity y of an ofF-shell gluon ex- 
changed in the f-channel can be written as 

^ /(qa, qb, y) = j d^-'q K(q„, q) /(q, q^, y), (2.1) 

where /(qa,q&,2/) is the gluon Green's function, and we have used boldface to denote 
transverse vectors. The BFKL kernel ^^(qi, q2) is split into contributions from virtual and 
real corrections, embedded in the trajectory a;(qf) and real emission kernel ifr(qi,q2), 
respectively 

i^(qi, q2) = 2 c^(q2) ^^-^(qi - qs) + i^,(qi, qs). (2.2) 

The first two terms in the perturbative expansion of each of these terms are known 

a;(q2) = a;(i)(q2) + u;(2)(q2) + . . . , KMi,^2) = i^(')(qi, q2) + i^(')(qi, q2) + • • • • 

(2.3) 

The leading contributions are given by 



with 



^ gl iVc r(l - e) 



The next-to-leading logarithmic correction to the trajectory is [1] 



27 ' 21N^ 



(2.6) 
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and to the real emission kernel [1] 
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(2.7) 



where 

L{z)=J^—Hl-t), C(n) = 5^fe-^ /5o=(^yiVe--n/J. (2.8) 

The contribution to the real emission kernel at NLL accuracy has three sources: 1) 
the one-loop corrections to one gluon emission in multi Regge kinematics, 2) two-gluon 
emission in quasi multi Regge kinematics, and 3) quark-anti-quark production in quasi 
multi Regge kinematics: 

i^(2)(qi,q2)=K(2).^/(qi,q2) + i^(2).fl9(q^,q2) + K(2).5%,,q2) (2.9) 

It is the purpose of this paper to calculate the contribution i^r^^''^^(qi, q2) in a form that 
retains not only the dependence on the sum of the transverse momenta of the emitted 
quark pair (qi — q2), but the dependence on the full momenta of both emitted particles. 

The divergences of i^r^''''^^(qi, q2) will fall into two categories; one at = (qi — q2)2 — 
regularised by the quark contribution to the trajectory at NLL, and an explicit pole for 
e — ^ 0, regularised by the quark contribution to iir^^^'^(qi, q2). The one-giuon production 
vertex for the collision of two Reggeiscd gluons at NLL accuracy can be calculated from a 
tree approximation for the quasi-multi- Regge kinematics combined with t-channel unitarity 
relations [4-6], or extracted directly from the one loop 5 gluon amplitude in the helicity 
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Figure 1: The quark-anti-quark production amplitude in Reggeon-Reggeon collisions (LHS) can 
be obtained as the sum of the contributing diagrams (RHS) according to the effective Feynman 
rules of Ref. [9]. The zig-zag lines denote the Reggeons (off-shell gluons). The BFKL evolution 
consists of a string of such vertices (and their LL and NLL gluon emission counterparts) connected 
with Reggeised gluon propagators. 



basis [7] . Expanded in I? = 4 -|- 2e it is given by [8] : 
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The divergent piece from the quark contribution is simply 
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(2.11) 
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2.1 The Quark-antiquark Contribution to the NLL Corrections 

The Reggeon-Reggeon-gg-vertex can be calculated using the Feynman rules of the effective 
action in the high-energy limit in Ref. [9]. The contributing diagrams are shown in Fig. 1. 
The benefit of this framework of obtaining the scattering amplitudes is that it gives an 
immediate picture of the contributing sub-processes. However, the RRqq (or g*g*qq) vertex 
can also be extracted directly from the full tree-level gg — > ggqq scattering amplitude in 
the helicity basis [10]. In terms of the operators of the effective action of Ref. [9], the result 
for the g*g*qq vertex jf^f^{qi,q2, h, ^2) is 



&Qi,Q2} = -gMki)[fHH{kuk2) - &f'b{k2,h)]v{k2) 



(2.12) 



where t*^ , are the colour group generators for the fundamental representation. In terms 
of the helicity amplitude A'^_^_ the spin, colour and flavour sum is 

E \lSfM,<l2,k,,k2)\' 

ii,i2,f 

2 gHN^ 
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(2.13) 
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and in terms of the light-cone momenta and complex transverse momenta the amplitude 
is given by [10] 
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where §12 = {ki + ^2)^ and £ = (^2 + 52)^ = (gi - ^1)' 



The quark contribution to the NLL corrections to the BFKL kernel, ifr^^'^^(qi, q2) in 
Eq. (2.17), is in the traditional analysis [1,2] obtained as the integral over the phase space 
of the emitted quark-antiquark pair of the fully exclusive RRqq-vertex 



1 



qi,q2j 



2qfqiiV2 
where k = {qi — q2 



dn dpf 6^^\qi -q2-ki- /C2) ^ \l'^^%{qi,q2, ki, k2] 

h,i2,f 

(2.15) 



^ = (^1 + ^2)^ = •S12 is the invariant mass, and 
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The result of the integration in D = 4 + 2e dimensions can be written as [11] 
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(2.17) 



This result enters directly in the real emission kernel of Eq. (2.7), and it is the validity 
of this result we will discuss. The problem we wish to highlight is the role of the S- 
functional in the integrations of Eq. (2.15). While each amplitude (e.g. 7j^^(q'i, ^2, ^1, ^2) 
and the gluonic counterparts) conserves energy and momentum (i.e. all 4 components of 

the momentum vectors), the conservation of the longitudinal and energy component is lost 
in Eqs. (2. 15), (2. 17), in order to write the result in a form that depends only on transverse 
vectors. In Eq. (2.17) there is no relation between the longitudinal and energy component of 



-5- 



the reggeon momenta, and the momenta of the quark-anti-quark pair. This has two effects. 
First, energy and momentum conservation is broken even at NLL accuracy when connecting 
the vertices to form the BFKL evolution (this effect was already highlighted in Ref. [3]). 
This part of the problem can be solved by using the framework of the direct solution of 
the BFKL evolution [3,12,13]. Secondly, within each Lipatov vertex the NLL corrections 
to the reggeon-reggeon scattering are not constrained to the correct energy probed in the 
corresponding LL process, but rather integrated over all energies. This obviously not only 
violates energy and momentum conservation at the level of each vertex, but also exaggerates 
the contribution of the NLL corrections. This problem, and its solution, is the focus of the 
remaining study. 

We will calculate the quark contribution to the NLL corrections by performing the 
integral in Eq. (2.15), but regularise the result by the phase space slicing method, and 
combine it with the pole in e of Eq. (2.11) from the quark contribution to the loop correc- 
tions of the one-gluon production Lipatov vertex. The result will be finite for ^ 0. In 
performing the integration, we will parametrise the phase space for fixed A in terms of the 
transverse momenta of the anti-quark (A;ix) and half the rapidity separation between the 
quark and the antiquark, r) (the quark and anti-quark are produced with rapidity r) and 
—77). We find 

1 dK dp/ (5(^)(gi -q,-k^-k2)= JdvJ (2.18) 

which once again illustrates that the NLL corrections to the Lipatov vertex includes contri- 
bution from a range of energies and not just that probed by the LL vertex (i.e. for a given 
qi , q2 the NLL corrections include contributions from all invariant masses and energies in 
the Reggeon- Reggeon collision, not just the one probed at LL). The parametrisation of 
the phase space in Eq. (2.18) allows for an immediate implementation of the integration in 
Eq. (2.15). 

Symmetry properties in the divergent part of the amplitude ensures that the integral 
over the l/N^ suppressed contribution of the quark-anti-quark vertex is finite, and the 
result for 



^(qi,q2) = 

64 (27r)^ / / 



(2.19) 



d^ki 
(27r)^ 



2 Re A^^_{h,k2)Af^{k2,h) - (|Af_(fci, A;2)|' + \Af_{k2,h)\'') 



can be directly compared to Eq. (23) of Ref. [11], which up to constants form the 
suppressed part of Eq. (2.17). The result is shown in Fig. 2 for fixed q2 = (0,15) GeV 
as a function of qi = (x, 0) GeV. We find complete agreement, after correcting a misprint 
in Eq. (23) of Ref. [11] (the denominator 216 qf should read 16 qf q^). The complete 
agreement verifies that our implementation of the amplitudes and parametrisation of phase 
space is correct. 
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Figure 2: The quantity J(qi,q2) of Eq. (2.19) for q2 = (0, 15) GeV and qi varying along (1,0) 
compared between the results obtained in this paper and the corrected results of Ref. [11]. 



2.2 Regularisation 

We now want to regularise the integral in Eq. (2.15) in order to implement the quark 
contribution to the fully exclusive Lipatov vertex at NLL. We will start by rewriting [10,14] 
the amplitude AV_{ki, 
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The quark-anti-quark contribution has two divergences: A singularity at k'^ — xA^ 
and one at A^ = 0. When constructing the regularised vertex, the first type of singularity is 
regularised by the corresponding singularity from the one-loop corrections to the one-gluon 
production, whereas the second type of singularity is regularised by the corresponding one 
in the gluon trajectory. In the framework of the direct solution to the BFKL evolution, the 
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regularisation of the latter divergence was performed already in Ref. [12,13]. In the present 
context we therefore need only to consider 7^ 0. In order to implement the exclusive 
NLL Lipatov vertex, we would need to perform the regularisation of the divergence at 
ki± — xA± using a phase space slicing method. Technically, we will perform the phase 
space integration by Monte Carlo techniques. If A'^_{ki, is probed for ki± — xA^ > A 
with A small (our reported results are stable for variations of A between 10~^ GeV and 
10~^*^ GeV), the standard result of Eq. (2.14) is used. When ki±_ — xA^ < A, in addition 
to the non-divergent pieces we return the average value of the integral of \A'^^_ ^jy(fei, ^2)!^ 
combined with the appropriate term from the one-loop quark contribution to the one-gluon 
emission vertex, evaluated for e — 0. 

It turns out that only the square of the two terms in Eq. (2.20) that have an explicit 
divergence at ki± x/S.± lead to a pole in e. All cross terms integrate to zero when 
integrated for feij^ — xAj^ < A. We therefore find 
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(2.22) 
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and the substitution ki — xA — 1 together with the following result 
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with L> = 4 + 2e, we can obtain the Laurent series in e for the relevant phase space integral 
of each of the contributing terms. The result is 



2\e 



A2 (27r) 
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3 r(l + e) £ 



(2.25) 



The generated divergence cancels with the divergent part of the one-loop quark contribu- 
tion to the one-gluon emission vertex in Eq. (2.11). 



3. Discussion 

Having thus constructed the regularised quark contribution, we can calculate the contribu- 
tion to the BFKL kernel as in Eq. (2.15) (remembering that also the divergent contribution 
from the virtual corrections are taken into account). The Monte Carlo method allow us 
to study directly the differential distribution in any variable. Of particular interest is the 
distribution in the energy of the Reggeon-Reggeon collision. The LL contribution come 
from the emission of a single gluon of energy |qi — q2|- However, when the NLL corrections 
are calculated as in Eq. (2.15), contributions arise from all energies above this scale. In 



-8- 




90 100 
E [GeV] 



Figure 3: The energy distribution of the contribution to the NLL corrections to a Reggeon-Reggeon 
collision of \/2 • 20 GcV in the standard approach. 



Fig. 3 we have plotted the distribution in the energy of the Reggeon-Reggeon collisions 
contributing to the NLL corrections to the LL result. We have chosen qi and q2 to be 
perpendicular, with |qi| = |q2| = 20 GeV, however the findings have general validity. The 
average energy of the Reggeon-Reggeon collisions included in the NLL corrections calcu- 
lated according to Eq. (2.15) is ~ 40% larger than the energy in the LL setup (for all angles 
and for an arbitrary value of |qi| = |q2|)- Furthermore, the average rapidity separation 
between the quark and the anti-quark is « .56 units of rapidity. Often, the contribution 
supposed to described the radiative corrections to a single gluon emission would be de- 
scribed as a two-jet configuration! All the quark-anti-quark configurations included in the 
integrations of Eq. (2.15) must be included in the evolution to obtain overall NLL accuracy. 
But assigning them all to a single transverse scale will potentially greatly exaggerate the 
impact of these NLL corrections on the evolution^. 

However, this is a problem that cannot be resolved as long as only transverse degrees of 
freedom are considered in the governing BFKL equation. We therefore propose a framework 
which will perform the BFKL evolution fully differentially in the momenta of all emitted 
particles, i.e. not only using the fully differential BFKL equation [3, 12, 13] 



-TT- ^7— /(qa,{fci±,yi},q6,y) = K qa + Vfci±,qa + VA;i± /(qa,{fci±,yi},qfe,y), 

(3.1) 



A similar problem will appear in the calculation of the impact factors 



but also keeping the NLL corrections fully exclusive as detailed in this study for the quark 
contributions. This will ensure not only the conservation of all components of the mo- 
mentum vector, but also that for a given set of Reggeon momenta, only the relevant NLL 
corrections are included. The result of such evolution would be in complete agreement with 
the description of the Reggeisation of t-channel gluon amplitudes according to Ref. [15], 
when also here the NLL vertices are kept exclusive. Of course the result of the standard 
BFKL formalism can be obtained by suitable phase space integrations and ignorance of 
the conservation of longitudinal momentum. 

The large spread in energies of the contributions to the NLL corrections will have 
sizeable effects in a realistic application of the BFKL framework. In the application of 
the BFKL framework to the description of the small- a:; behaviour of partons, the calculated 
NLL corrections to the kernel cannot be ascribed to a single value of x. And in the emerging 
framework of applying the BFKL evolution in the description of the production of multiple 
hard jets, the NLL corrections would be weighted with different pdf-factors, according to 
their relevant contribution to the centre of mass energy. This effect will change the impact 
of the NLL corrections, and could further reduce their importance. 

4. Conclusions 

We have presented a new calculation of the quark contribution to the NLL corrections to 
the BFKL kernel, based on an explicit Monte Carlo integration of the produced particles. 
We have demonstrated that in the standard calculation of the NLL BFKL kernel, the NLL 
corrections arise from Reggeon-Reggeon energies far higher than that which is probed at 
LL. This exaggerates the apparent effect of the NLL corrections. 

We have presented a method which will assign the NLL corrections to the appropriate 
LL contributions. We are looking forward to reporting on the gluon contribution to the 
NLL corrections, and on the implementation of the new framework in the prediction for 
multi-jet events at the LHC. 
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